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, The influence of A isobar components on the ground state properties of nu- 

q^I clear systems is investigated for nuclear matter as well as finite nuclei. Many- 

body wave functions, including isobar configurations, and binding energies are 
evaluated employing the framework of the coupled-cluster theory. It is demon- 
strated that the effect of isobar configurations depends in a rather sensitive 
way on the model used for the baryon-baryon interaction. As examples for 
realistic baryon-baryon interactions with explicit inclusion of isobar channels 
\Q \ we use the local (V28) and non-local meson exchange potentials (Bonn2ooo) 

but also a model recently developed by the Salamanca group, which is based 
on a quark picture. The differences obtained for the nuclear observables are 
related to the treatment of the interaction, the 7r-exchange contributions in 

| particular, at high momentum transfers. 
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I. INTRODUCTION 

> : 

The development of efficient computing facilities has enabled very sophisticated calculations 
for the solution of the nuclear many-body problem. Starting from realistic models for the 



nucleon-nucleon (NN) interaction, which give very accurate fits of the NN scattering data 
below the threshold for pion production JI|-§[, one can solve the few-nucleon problem up 
to A = 8 nucleons in a way which yields essentially the exact solution Q. Introducing an 
additional three-nucleon force |||6| one can obtain results for the basic low-energy properties 
of these nuclei, which are in good agreement with the experimental data. 
This demonstrates that the low-energy properties of nuclei are well described within the 
conventional model of nuclear physics, in which nuclei are considered as a system of nu- 
cleons, treated as inert particles interacting via two-body forces. All sub-nucleonic degrees 
of freedom, which may lead to modifications of the hadrons in the nuclear medium, and 
dynamical relativistic j7| effects are represented by a phenomenological three-nucleon force. 



*On leave from University of Salamanca, Spain. 
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On the other hand, however, one knows that nucleons cannot really be considered as el- 
ementary particles and sub-nucleonic degrees of freedom, like e.g. the possibility to excite 
strongly interacting nucleons, could be very important. In particular the excitation of nu- 
cleons to the A(3, 3) resonance may have some effect on the low-energy and bulk properties 
of nuclear systems. First investigations on the importance of isobar degrees of freedom have 
been performed more than 20 years ago |8|-p|. Those studies demonstrated that isobar 
configurations yield an important contribution to the medium range attraction of the NN 
interaction. Conventional models of the NN interaction account for these mutual polari- 
sation of the interacting nucleons in a phenomenological way. For example, a part of the 
a meson exchange in One-Boson-Exchange (OBE) models for the NN interaction can be 



related to such isobar terms 12 



In a conventional nuclear structure calculation this part of the NN interaction is identical in 
the nuclear medium as compared to the vacuum where the effective NN interaction has been 
adjusted to describe the NN scattering data. If, however, the isobar degrees of freedom are 
taken into account explicitely, one obtains a modification of the NA and AA propagator in 
the medium. This implies that the effective NN interaction including such intermediate iso- 
bar states is different in nuclear matter as compared to the vacuum case. The corresponding 
part of the medium range attraction is quenched. This feature has been investigated by var- 
ious groups using the Brueckner-Hartree-Fock approximation §-i [| or within a lowest order 



variational calculation |L3| and binding energies were obtained, which were much weaker 
than in corresponding calculations ignoring the explicit treatment of isobar excitations. 
For nuclear matter at higher densities the explicit consideration of isobar configurations leads 
to an enhancement of the pion propagator which has been called a precursor phenomenon to 



a phase transition of pion condensation [14-16]. This leads to rather attractive contributions 



to the binding energy which originate from ring diagrams involving A-hole excitations [|I7| . 
Nuclear structure studies including isobar excitations have furthermore been performed for 
few nucleon systems |T8|-|2~ 



Most of these older studies have been performed using rather simple models for the baryon- 
baryon interaction. The transition potentials describing the NN — > NA and NN — > AA 
were approximated in terms of local 7r-exchange potentials. During the last years new models 
for the baryon-baryon interaction have been developed. It has been demonstrated that a local 
approximation of the 7r-exchange term tends to overestimate these contributions considerably 
pl| , |22[| . This is also true for the transition potentials leading to isobar excitations ||23|| . So 



it is one aim of the present investigation to update nuclear structure studies with explicite 
treatment of isobar excitations using modern models for the NN interaction. We are going 
to compare results for nuclear matter and finite nuclei, calculated for the Argonne V28 
potential p4]| , a recent update of the non-local meson-exchange potential denoted as Bonn 20 oo 
|25fl and a model which has recently been developed by the Salamanca group [p6p . This 



Salamanca interaction is derived in the framework of the Chiral Quark Cluster (CQC) 
model. The problem of two interacting clusters (baryons) of quarks is solved by means 
of the resonating group method. The Pauli principle between the interacting quarks is an 
important source for the short-range repulsion of the NN interaction | 27fl . At large distances 
the 7r exchange between the quarks in the two clusters evolves to the n exchange between 
two baryons. At shorter distances, however, this non-local model for the baryon-baryon 
interaction might yield results that are quite different from a meson-exchange picture. This 
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Salamanca potential does not give such a perfect fit to the NN scattering phase shifts as 
the Bonn 2 ooo or the V28. For the 1 Sq, 1 Po and 3 Si — 3 D\ partial waves of the NN system, 
however, the agreement with the empirical data is rather good. 

An explicit evaluation of isobar components in the nuclear wave function is also motivated 



from recent experiments, which try to measure such isobar components |28 . 
The isobar components in the nuclear wave function and the resulting ground-state proper- 
ties will be evaluated in an extension of the coupled cluster method []29| ]. This extension is 
presented in section 2 where we will also compare predictions of the coupled-cluster method 
with the Brueckner-Hartree-Fock approximation. Results for the binding energy and isobar 
probabilities obtained for nuclear matter and finite nuclei will be presented in section 3. 
Special attention will be paid to the difference between the various interaction models. 



II. COUPLED-CLUSTER APPROACH WITH ISOBAR EXCITATIONS 

In the coupled-cluster approach |2P| one starts assuming an appropriate Slater determinant 
$ as a reference state for the system under consideration. In the examples considered below 
this reference state will be a Slater determinant defined in terms of appropriate oscillator 
single-particle wave functions for the case of 16 O, while for the case of infinite nuclear matter 
$ stands for the antisymmetrized wave function build in terms of plane waves with momenta 
less than the Fermi momentum kp. The exact eigenstate \I/ is then written as 

^ = e 5 $, (1) 

with 5* being an operator of the form 

S=J2 Sn, (2) 

71=1 

where S n is an n-particle operator and in order to be complete one has to consider operators 
up to n = A with A the number of baryons in the system. The operator S n describes the 
formation of an n-particle n-hole excitation relative to the reference state $. For the case 
of n = 2 it can be written 

S 2 = - (PiP2|5 , 2|t / i^2> a ] Pl a ] P2 a U2 a Ul . (3) 

In this equation a p . stand for fermion creation operators in states which are unoccupied in 
$, while a Vi represent annihilation operators for the nucleon single-particle states which are 
occupied in the Slater determinant <3>. Note that the aj,. may also represent the creation 
of A isobar states. Therefore the S 2 amplitudes describe two-particle two-hole excitations 
relative to $ but also NA and AA excitations. 
One can now use the Schrodinger equation in the form 

e- s He s <$> = E<$> , (4) 

and project this equation on the reference state $ and n-particle n-hole states relative to $ 
which we will identify by & Pl ...p nVl .. Mn - This leads to an expression for the energy 
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E = ($\e- s He s \$ 



(l + Si + ^ + S^ \$\ ) ( 5 ) 

and to a set of coupled equations for the amplitudes of linked n-particle n-hole excitations 
S n . This set of equations is truncated by assuming that amplitudes S n with n larger than a 
given value m can be ignored. As an example we consider the m = 2 approximation, i.e. we 
ignore the effects of linked 3-particle 3-hole excitations and higher, and in order to simplify 
the notation, we furthermore assume that we have chosen the reference state such that Si 
vanishes (note that this is true in particular for infinite nuclear matter because of the the 
translational symmetry). In this case we can write the correlated two-body state as 

>A= (1 + S 2 )\l>iV 2 >A , 

where a subscript A is used to identify antisymmetrized two-body states. With these sim- 
plifications, the equation for the amplitudes S2 can be reduced to 

(pip 2 \{Ti + T 2 )S 2 \viV 2 ) A ~ Yl {(PiP2\S 2 \viS2) A (v\h\vi) + (pip 2 \S 2 \viv) A (v\h\v 2 )} 

V 

+ (pip2\Vx2\viV2) A + (P1P2IS2I vv') A (vv'\Vx 2 \viV 2 ) A = 



(*) 

(6) 

Note that the term identified with (*) includes a summation over intermediate hole states. 
If we ignore this term and use furthermore that the single-particle hamiltonian of nuclear 
matter is diagonal in the plane wave states, (a\h\(3) = e a 5 a p, we can rewrite @ into 

Vx2Wi^) A = V\viv 2 ) A + V ■ ^= 7^X2^2) A 

e vi t € U2 — ii — i 2 

S v ' 

= V\uiu 2 ) A + VS 2 \viv 2 ) A . (7) 

In this equation we have introduced the starting energy u and the Pauli Operator Q P , 
which restricts the sum over intermediate states to those of the form \pip 2 >. This means to 
nucleon single-particle states which are unoccupied in the reference state $ or to isobar 
excitations. If we identify V\2 with the Brueckner G-matrix, the equation (|7|) takes the 
form of the Bethe-Goldstone equation 

G (u ) = V + V-%-G(u) , (8) 

with H the operator of the kinetic energy, i.e. assuming the conventional choice for the 



spectrum of the particle states |22|j3(]|j3T 



In order to visualize the relevance of the hole-hole term (*) in (H) we have performed cal- 
culations for nuclear matter with and without this term, using two different models for the 
NN interaction, which do not include isobar degrees of freedom, explicitely. One of these 
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examples is the Argonne V14 potential f24]| , which is defined in terms of 14 operators, each 
of them multiplied with a local potential. The second example is the charge dependent 
Bonn potential |]J, a meson exchange interaction, which is evaluated in momentum space 
and contains non-local contributions. 

The results of the energy of nuclear matter as a function of the Fermi momentum hp, 
calculated in the Brueckner-Hartree-Fock (BHF) approximation is displayed by the dashed 
lines in Fig. III. The differences originating from the two interaction models have been 



discussed e.g. in [^TJ: A local interaction (V14) tends to be stiffer than a NN interaction 
based on the non-local meson exchange model (CDBONN), fitting the same NN scattering 
data. Since the two-body correlations in nuclear matter are quenched as compared to the 
case of NN scattering in the vacuum, a stiffer interaction tends to predict less binding energy 
than a softer one. As a consequence, the BHF energy calculated for the V14 interaction is 
much less attractive than for CDBONN. 

If also one includes the hole-hole ladder terms in solving eq.(Q), one obtains the corresponding 
solid lines in Fig. |l|. The comparison shows that the effect of these hole-hole ladders on the 
calculated energy of nuclear matter is rather weak in this range of densities. The hole-hole 
ladders yield an effect which is weakly repulsive. The effect is a little bit larger for the V14 
interaction as compared to the softer CDBONN. All these results indicate that the coupled 
cluster approach, restricted to the S m amplitudes with m < 2 yields results very similar 
to the BHF approximation, employing the conventional choice for the intermediate particle 
spectrum. 

However, it is not the issue of the present work to perform nuclear structure studies within 
the conventional approach. More sophisticated calculations including up to three-hole line 
terms in the Brueckner expansion scheme for nuclear matter [31] or coupled-cluster calcula- 
tions for finite nuclei including S3 terms have been performed P2"| , |3"3"|] . 
The central aim of this work is to account for the isobar excitations using the coupled- 
cluster approach restricted to m less or equal 2. For that purpose we consider the baryon 
baryon interaction models V28 P3f| , Bonn 2 ooo PH and the Chiral Quark Cluster (CQC) 



model developed in Salamanca ||26|| , which all include the scattering to NA and A A states. 
While V28 and Bonn 20 oo yield rather accurate fits of the NN phase shifts in all partial waves, 
the CQC model leads to such a good fit only for the channels with isospin T=0, J=l and 
T=l, J=0 channels. Therefore we have replaced the CQC model by the Bonn2ooo interaction 
model in all other channels. 

For the case of nuclear matter the equations (0) has been solved as an integral equation 



employing the usual angle-average approximation for the Pauli operator [gj] . The hole- hole 



ladder term can be introduced as an additional non-linear term, for which self-consistency 
is obtained in an iterative procedure. 

In the case of finite nuclei we have solved the coupled- cluster equation by considering an 
expansion of the correlated two-body wave function in a basis of relative wave function 
defined in a box of a given radius. This basis provides an independent control of the max- 
imal distance and momentum relevant for the correlated waves. The method, restricted to 



nucleon-nucleon correlations only has been described in ||35|| . The extension to include isobar 
configurations is straight forward p3[ . 

As a result of the calculation we not only obtain the energy of the system ([5]) but one can 
also calculate the one-body density matrix 
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using the techniques described in [[36]. This one-body density allows for the evaluation of 
the radius of finite nuclei and also yields the probability that a A is excited. 



III. RESULTS AND DISCUSSION 



The energy per nucleon calculated for homogenous nuclear matter at various densities is 
displayed in Fig. [|. The results obtained for the three different interactions treating isobar 
excitations explicitely (Argonne V28, Bonn 2 ooo an d Salamanca CQC) are compared to those 
obtained within the conventional framework using the Argonne V14 interaction model. All 
results have been obtained in the coupled-cluster or exponential S approach restricting the 
excitation operator to S m with m < 2. All calculations including isobar configurations yield 
results for the binding energy, which are less attractive than the result obtained for the 
conventional calculation. 

The reason for this loss of binding energy has been presented already long time ago (see 
e.g. |H|) and we just want to repeat it using the language of perturbation theory. The 
contribution of isobar configurations to the effective interaction of 2 nucleons can be written 
in lowest perturbation theory as 

= 4a %— V na + Vl A -— V^a , (10) 

where Vna and Vaa represent the transition potentials for the NN — > NA and NN — > AA 
transitions, respectively. The energy of the interacting nucleons is denoted by uj and H NA ( 
//aa describe the hamiltonian for the intermediate NA (AA) states. These contributions 
are responsible for a sizeable part of the medium range attraction of the NN interaction. 
Therefore a realistic interaction model like the Argonne V14, which does not allow for isobar 
configurations contains attractive components which simulate the effects of isobar excitations 
in fitting the NN scattering phase shifts. The effects of the attractive isobar terms displayed 
in (|1(]) are reduced in the nuclear medium, since the interacting nucleons are bound (u 
gets negative) and because a part of the NA configurations is unavailable due to the Pauli 
blocking. This quenching of the attractive isobar terms is observed only, if the isobar effects 
are treated explicitely, it is not contained in the conventional models, that simulate the 
isobar effects in the effective NN interaction in a pure phenomenological way. 
These arguments explain the loss of attraction due to the explicit inclusion of isobar ex- 
citations. Since the reduction of the isobar terms increases with density, they can also 
explain that the repulsion increases with increasing density, a feature which tends to shift 
the saturation point to lower densities. 

The calculated binding energy is smallest for the Argonne V28 interaction model and slightly 
larger for the Bonn 20 oo and the Salamanca CQC model. To some extent this could be ex- 
plained by the observation that local interaction models, like Argonne V28, tend to predict 
a weaker binding than non-local interactions, which fit the same NN phase shifts ||21|| . This 



feature, however, may also be interpreted as an indication that the predicted isobar effects 
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are larger for Argonne V28 than for the other two models under consideration. This inter- 
pretation is supported by the calculated probabilities of isobar excitations in nuclear matter, 
displayed in Fig. |[ 

The prediction for the isobar probabilities derived from the various interaction models differ 
in a very significant way. At the empirical saturation density, which corresponds to a Fermi 
momentum kp = 1.36 fm _1 , the difference is larger than a factor two (see also Table |). De- 
spite these differences in the prediction of the total A probability, there are some common 
features in the predictions of these models with quite different origin. If one tries to analyze, 
which partial waves provide the most important contributions, one observes that all inter- 
action models predict a larger contribution from the excitation of AA configurations, which 
can occur in interacting pairs of baryons with T = and T = 1, than from the excitation of 
iVA excitations, which occur in T = 1 partial waves only. 

All interaction models predict large contributions from those partial waves, in which the 
interacting nucleons are in a state with relative angular momentum I = 0. However, for a 
complete calculation one cannot ignore the contributions from the higher partial waves (see 
also Table |). 

The features we have discussed so far for the case of infinite nuclear matter are also observed 
in the results obtained for the finite nucleus 16 0, which are displayed in Table |TJ. The 
calculated binding energy per nucleon is smallest for the Argonne V28, about 0.7 MeV 
and 2.4 MeV per nucleon larger for the Bonn 20 oo and Salamanca model, respectively. The 
increase in the calculated binding energy is correlated with a smaller prediction for the A 
probability Pa- Also for 16 O we observe a difference by almost a factor of 2 between Pa 
derived from V28 and Salamanca CQC model. The probabilities Pa calculated for 16 O are 
similar in magnitude as those evaluated for nuclear matter at small densities (kp around 1 
fm _1 ). Therefore a nuclear matter calculation of isobar effects seems to provide a reasonable 
first estimate for the case of finite nuclei, if one uses a local density approximation. For a 
more detailed information, like the relative importance of different partial waves, an explicit 
calculation of the finite systems is required. 

The differences in the calculated energies and A probabilities obtained in these different 
interaction models originate of course mainly from the different transition potentials V/va 
and Vaa (see also eq.fllQf)). Some of these differences, like the treatment of the tt exchange 



contribution or the short-range behaviour in the local V28 approximation as compared to 
the non-local calculation in the other models have been discussed already in [p3|l . These 
differences are also the main origin for the different predictions obtained in the present 
calculation. As a typical example we show the amplitudes 



{NA 5 D \S 2 \Os 1/2 Os 1/2 ) j=ot=i , (11) 

calculated for 16 as a function of the relative distance r of the NA pair in the upper 
part of Fig. [|. The lower part contains the corresponding correlation function for the AA 
configuration. Inspecting the differences obtained from the three interaction models, one can 
see that the V28 interaction leads to a larger amplitude at larger distances r than the other 
two. This can be related to the fact that V28 uses a local n exchange contribution in the 
transition potential, which does not account for retardation effects which are due to the iVA 
mass difference. The different behaviour in the correlation functions S 2 at small distances r 
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must be related to the different kinds of models. While the short-range behaviour in the V28 
and Bonn 20 oo model are controlled by local and non-local form-factors, respectively, it is the 
coupled channel calculation within the chiral quark model, which provides the short-range 
behaviour in the Salamanca model. 



IV. CONCLUSIONS 

In the present investigation we try to compare the predictions for the bulk properties of 
nuclei derived from three different baryon-baryon interaction models, which account for 
isobar configurations explicitely. The main differences can be related to the models for 
transition potentials describing NN — > iVA and NN — > AA transitions. The quark model 
of the Salamanca CQC approach predicts weaker transition amplitudes at short range than 
the more phenomenological cut-offs in the Bonn 20 oo and Argonne V28 interaction. The long- 
range components of these amplitudes are dominated by the 7r-exchange, which is weaker in 
the non-local models (Salamanca and Bomi2ooo) than in the local interaction model (V28). 
These differences in the interactions are responsible for the differences in the predicted A 
probabilities in 16 O, which vary between 1.98 percent derived for the Salamanca model and 
3.87 percent for Argonne V28. The results are in fair agreement with the estimates derived 
from experiment [E8]. A large probability for isobar excitations is related to weak binding 



energy. Similar results are obtained for infinite nuclear matter. The isobar effects discussed 
here would correspond to the inclusion of a repulsive three-nucleon force in conventional 
nuclear structure calculations. 

Comparing the calculated binding energies with empirical values, one must keep in mind 
that the coupled-cluster approximation (including terms up to S 2 ), essentially corresponds 
to the Brueckner-Hartree-Fock approximation with a conventional choice of the intermediate 
particle spectrum. More binding energy is obtained from including three-body terms or us- 
ing the so-called continuous choice spectrum. This is known from corresponding calculation 
using nucleon degrees of freedom only [[31 1. The explicit treatment of isobar configuration, 



however, would give rise to additional three-body terms, which are not yet taken into ac- 
count. The effects become very large at higher densities if one considers old local models for 



the isobar excitations [pT6| JlTfl . The isobar effects may provide reasonable corrections if the 
modern interaction models are considered. 
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v 28 


Bonn 2 ooo 


CQC 


P A (total) 


8.67 


6.55 


4.00 


P(NA) 


4.03 


2.79 


1.87 


P(AA) 


4.64 


3.76 


2.13 


Pa(%) 


2.77 


2.43 


1.24 


Pa( 3 Si) 


1.80 


2.00 


0.64 



TABLE I. A probability per nucleon in nuclear matter at saturation density derived from 
various interaction models. The total probability originates from the excitation of iVA (P (NA)) 
and AA configurations. Also shown are predictions if only the coupling to A configuration from 
NN channels 1 Sq and 3 Si are considered. All entries are in percent. 





v 28 


Bonn 20 oo 


CQC 


E/A [MeV] 


-2.73 


-3.49 


-5.15 


r [fm] 


2.81 


2.65 


2.53 


e sl/2 [MeV] 


-34.65 


-37.88 


-44.25 


£ P 3/2 [MeV] 


-16.49 


-18.54 


-22.03 


e pl/2 [MeV] 


-13.56 


-14.82 


-17.67 


Pa [%] 


3.87 


3.71 


1.97 


Pa (NA) [%] 


1.75 


1.43 


0.88 


Pa (AA) [%] 


2.12 


2.28 


1.08 



TABLE II. Energy per nucleon (E/A), radius (r), single-particle energies for the nucleons and 
A probability per nucleon for 16 O. The total probability originates from the excitation of NA 
(P(iVA)) and AA configurations. Results are presented for the various interaction models dis- 
cussed in text. 
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-20 1 ' 1 ' 1 ' 1 ' 1 

1 1.25 1.5 1.75 2 

Fermi momentum [fm '] 

FIG. 1. Binding energy of nuclear matter as a function of the Fermi momentum. Results are 
given for the Argonne V14 and the CDBONN interaction, using the BHF approximation (dashed 
lines) and with inclusion of hole- hole ladder terms (solid lines) . 




FIG. 2. Binding energy of nuclear matter as a function of the Fermi momentum. The contri- 
butions from hole-hole ladders, which are negligibly small (see Fig. 1), have been ignored in the 
results displayed here. Results are given for the various interaction models with explicit considera- 
tion of isobar excitations. For a comparison the result obtained for the conventional Argonne V14 
interaction is also included. 
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1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

1 1.2 1.4 1.6 1.8 2 

Fermi momentum [fm ] 



FIG. 3. A probability per nucleoli in nuclear matter as a function of the Fermi momentum. 
Results are presented for the various interaction models discussed in text. 




12 3 4 



r [fm] 

FIG. 4. iVA (upper part) and AA correlation function originating from two nucleons in the 
0si/2 shell of 16 O as a function of the relative distance r. 
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